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Abstract We extend the Poincare-Lyapounov-Nekhoroshev theorem from 
torus actions and invariant tori to general (non-abelian) involutory systems of 
vector fields and general invariant manifolds. 

Introduction. 

The celebrated Poincare-Lyapounov theorem gives conditions ensuring that a 
periodic solution of a smooth dynamical system is persistent under small per- 
turbations. 

The theorem was extended by Nekhoroshev pp to the case of quasi-periodic 
solutions of partially integrable hamiltonian systems. His result is referred to 
as the Poincare-Lyapounov-Nekhoroshev (PLN) theorem. 

Detailed proofs of Nekhoroshev's result were provided in [S] from an ana- 
lytical point of view, and in |17j from a geometrical one (the latter work also 
contains an extension to non-hamiltonian vector fields). The non-hamiltonian 
frame was fully considered in JHjj where generalization of results concerning 
bifurcation from periodic solutions to bifurcation from quasiperiodic solutions 
are also considered. 

Further developments include extension to infinite dimensional systems (and 
existence of breathers) [7], to perturbation of systems with non-compact invari- 
ant manifolds |141 ITS] , and to partially integrable bi- hamiltonian systems [2U] ■ 

All these results are based on the assumption that the algebra Q of vector 
fields under considerations, hamiltonian or otherwise, is abelian; and correspond- 
ingly the invariant manifold A is a torus T k or (in the non-compact case 
the product of a torus by a contractible manifold. 

After the publication of |17j prof. Duistermaat remarked in a kind letter 
that, by some arguments based on the geometry of foliations, one should expect 
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an equivalent result to hold also in the non-abelian case. 1 The purpose of this 
note is precisely to extend the PLN theorem to the case of non-abelian algebras 
of vector fields, and more generally to involutory systems of vector fields, albeit 
in a slightly different way. 

Some words are maybe in order, before going into mathematical detail, about 
the physical motivation for such a study and possible physical applications of 
its results. 

The main field of applications for the results obtained here would be that 
of (differentiable) dynamical systems - i.e. systems of first order ODEs on a 
smooth manifold. In this framework, one of the considered vector fields would 
be the dynamical one, while the other ones would be symmetries of the former 2 . 
It should be recalled, indeed, that for first order (systems of) ODEs, contrary to 
all other cases, the natural algebraic structure for the set Q of vector fields being 
Lie symmetries of this is not that of a Lie algebra (Q being infinite dimensional 
as such EHj) but rather that of a Lie module, Q being finite dimensional 
as such ^] In geometrical terms, this corresponds to having a finitely 

generated set of vector fields in involution a la Frobenius. 

As detailed below, our results is of interest mainly when the invariant man- 
ifold whose persistence is considered has a nontrivial topology. It is well known 
- and rather obvious - that a symmetry of a vector fields maps solutions with 
a given topology into solutions with the same topology (we are here referring 
to the topology of trajectories for the solutions) m^I- The result we give 
here can also be seen as a generalization of this to the case where the vector 
fields as well as the symmetries depend on control parameters, and moreover to 
encompass also the case where the symmetry vector fields also have (at least a 
set of) trajectories with compact closure: this sets further restrictions on the 
persistence of the resulting compact invariant manifolds. 

Acknowledgements. This work was triggered by the remarks (a rather long 
time ago) of prof. J.J. Duistermaat on my previous work 17 . I also grate- 
fully acknowledge useful discussions with N.N. Nekhoroshev as well as with D. 
Bambusi, P. Morando and J. Pejsachowicz. 

1 Statement of the problem 

In this section we describe the general setting of our problem, i.e. persistence of 
a (/-invariant submanifold Ao where Q is an involutory system of vector fields. 

We will not use the most general setting, but the most physically relevant: 
we deal with a system of parameter-dependent vector fields on a given manifold 

x In the appendix of 1171 it was remarked that several parts of the proof of the PLN theorem 
given there do not extend to the case where Q is non abelian and A is not a torus. The result 
we obtain here is indeed weaker than the one holding in the abelian case, and the proof 
requires some modification of the arguments used there. 

2 Note that considering these on the same footing corresponds to what is done in the 
hamiltonian case, where the dynamical Hamiltonian and those describing the commuting 
integrals of motion are treated on equal basis. 
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V (the phase space), so that M is the direct product of the phase space and a 
parameter space Q (a discussion of the general case will be given elsewhere). 

Let V be a smooth manifold of dimension p, and Q = {Xi, Xj} a set 
of smooth vector fields in 7-\ which can depend on external parameters e, say 
e G R«. 

We assume (i) that for e — there is a smooth submanifold Ao C V which 
is invariant under Q, i.e. such that Xi : Aq — > TAo for all i = 1, d. 

We assume moreover that there is a tubular neighborhood V of A such that, 
for |e| small enough, (ii) the vector fields Q are in involution - in the Frobenius 
sense - in V, and (iii) span a regular distribution in V. 

Let us briefly recall what these assumptions mean. The Xi are in Frobenius 
involution if [Xi,Xj] = afj(x)Xk, with ofj G C°°(V A ) smooth real functions 
(depending also on the external parameters e) on V. 

Also, denote by V x C T X V the distribution associated to Q at the point 
ieP; that is, 

:= {veT x V : u = X\ x , X G £} . 

(Note that S C 7 C P is ^-invariant if, for all x G S, V x C T x 5.) Then 
£> = Up Z? x is the distribution on TP associated to Q\ this is regular in U if 
the subspaces V x C T X V have constant dimension for ieV. 

We would like to identify conditions ensuring that A is part of a smooth 
family of ^-invariant manifolds A e , isomorphic to Ao. 

It will be convenient to consider the product space Af = V x R 9 , and corre- 
spondingly U — V x Q with Q a suitably small neighborhood of the origin in 
R 9 . Our problem can be studied in U. We will also write N c for the intersection 
of M with the level manifold e = c; by construction, M c — A/o = "P. 

2 The fiber bundle construction 

We will see U as a fiber bundle (U, ir, Ao) over Ao, with contractible fiber. Our 
problem amounts then to the problem of identifying conditions which ensure 
the existence of a family of C/-invariant near-zero sections of this bundle. 

We are led to consider a certain (in general, nonlinear) connection V on this 
bundle, and covariantly constant sections of the bundle under V. 

We define for any point x G Ao a local smooth manifold E x C U of dimension 
s = n — k which is transversal to Ao, with Yj x n Ao = {x}, so that no two such 
manifolds intersect, and they define a smooth distribution in U 3 . We also 
consider linear local manifolds S x tangent to Y. x in x (the manifolds S x and 
S x are canonically identified by V); or to choose S x to be linear, which is fully 
legitimate PS] , 

In this way U is a bundle over Ao, and T, x represents the fiber through 
x G Ao; we denote the corresponding projection as ir, and write the bundle as 
(U, 7r,Ao) with 7r _1 (x) = T, x . 

3 If M is equipped with a riemannian metric, we can take as the local geodesic manifold 
through x orthogonal to Aq in x. 
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We can also consider Uo — U C\Nq — V, and a x := £ x n£/o; these are smooth 
manifolds of constant dimension s = n p — k (and codimension n q in ~S X ), and 
we denote by s x = S x D Uq the corresponding local linear manifolds. Then Uq 
is also a bundle over A , and a x represents the fiber through x G A ; we denote 
the corresponding projection as ttq, and write the bundle as (Uq,tto, Aq) with 
7Tq 1 (x) = a x . 

The distribution T> associated to Q has constant dimension r in U , by hy- 
pothesis. Moreover, V is tangent to Ao, and thus transversal to S a in x for 
all x G A; this implies that T> is also transversal to E x for all points u G T, x 
sufficiently near to i. If U has been chosen to be sufficiently small - which we 
assume from now on - the transversality condition is met for all points in U. 
Note that V defines a canonical identification between the local manifolds Y, x 
and S x defined above. 

We assumed moreover that k — r, hence V defines a (Frobenius integrable) 
distribution of horizontal spaces in U and thus a (in general, nonlinear) con- 
nection V in U EZ| • As this is defined by G, it will be referred to as the 
^-connection in U. 

The problem of existence of smooth (/-invariant manifolds isomorphic to Ao 
and near to it is, with this construction, translated into the problem of existence 
of near-zero (/-invariant sections for U, i.e. of sections of U which are invariant 
under the ^-connection V. 

Note that V-invariant sections always exist locally over any chart A in Ao, 
due to Frobenius theorem, since Q is an involutory system; thus the nontrivial 
part of the problem is purely global. 

As we deal with a small neighborhood of A we could - and we will indeed 
- consider the (transverse) linearization of the Xi around Ao. Correspondingly, 
we can consider the (transverse) linearization Vo of the connection V. 

We stress that V is in general a nonlinear connection, acting nonlinearly on 
sections (this is not a problem, as we are only interested in fixed points under 
this action), while Vo is - by definition - linear and has an associated covariant 
derivative acting linearly on sections. Finally, we note that, as it follows at once 
from Q being Frobenius integrable, the connections V and Vo in U are (locally) 
flat. 

3 Loops, Poincare-Nekhoroshev map, and mon- 
odromy matrices 

Let us fix a reference point x G Ao, and a loop 7 through x in Ao. Consider 
then a point w G 7r _1 (x) = S x ; the loop 7 is lifted by the (/-connection to a 
curve 7 in U, which defines a map in S x . 

We associate in this way a (monodromy) to each loop 7 

in Ao. These maps obviously form a group, the holonomy group TL X at x. When 
we fix a basis for the homology of A, we will refer to the maps {A^ 7l , -M 7r } 
as a set of generators for H x . 
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When we consider only contractible loops, the corresponding subgroup TL X C 
7i x is the restricted holonomy group at x. This is a normal subgroup in 7i x , 
and H x /H x is discrete [21] ■ It follows from our assumption that Q is Frobenius 
integrable (via the Ambrose-Singer theorem 2 24J) that the restricted holonomy 
group H x is trivial, H x — {e}. 

Now, for x an arbitrary reference point on Ao and E x = ir~ 1 x, V-invariant 
sections of [/ correspond to points w G S x which are fixed points for all elements 

Later on, we will find convenient to focus on the linearization of M 7 : T, x — > 
S K at a;; this is a linear operator acting in the linear space V = T X T, X ~ R n , 
i.e. a (n — fc)-dimcnsional matrix, called the monodromy matrix (at x) for 
the loop 7. We denote this as M T 

In this respect, recall that M. 1 also defines a map V 1 : S x — » S x . This will 
also be called the Poincare-Nekhoroshev map for 7 If we identify 

S x with a neighbourhood of the identity in V, as can be done via V, then M 7 
is exactly the linearization of this map. 

Monodromy matrices for all loops through x clearly form a group, called 
the monodromy group at x, and representing the linearization of the holonomy 
group at x. Thus we also refer to this as the linearized holonomy group, and 
denote it as h x . Given a set of generators for 7i x , the corresponding linearized 
operators will be a set of generators for h x . 

Remark 1. As well known, monodromy groups based at different points are 
conjugated, and homotopic loops based at the same point provide the same 
monodromy maps and matrices. 

4 The difference Poincare-Nekhoroshev map and 
its linearization 

We have introduced the Poincare-Nekhoroshev map "P 7 : S x — > S x , and its 
linearization around A, i.e. the monodromy matrix Af 7 :V—fV. By definition, 
V-yX = x; at the linear level, x corresponds to the origin in V, and M 7 Q = by 
linearity. Invariant near-zero sections are obtained for the near-zero £ £ V such 
that M 7 £ = £ for all loops 7. 

Instead of considering V~ t and M 7 , it is more convenient (as in the abelian 
case) to consider the map IZ-y := J — : S x — > S x and its linearization (also 
called the linearized difference Poincare-Nekhoroshev map) i? 7 = I — Af 7 ; note 
i? 7 : V -» V 4 

Let us now focus on a given x G Ao and a given loop 7 C Ao through x, and 
look for ^ £ Sj such that 7?. 7 £ = 0. By definition, £ = x satisfies this, and we 
are interested in knowing if there is any nearby £ satisfying this equation. 

4 With the construction considered in 1171 1181 we had to quotient out the action of Q from 
the Poincare-Nekhoroshev map; this is not needed now as we are directly considering the 
5-connection and not the flow under specific vector fields X 6 Q. Actually, M.^ could be 
seen as a map between local sections (invariant under the (7-connection) , see the discussion in 
(T7IIT8) . 
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This leads to investigate the question of existence of any zero eigenvalue for 
the linear map i?y if this is the case, the zero eigenspace corresponds to fixed 
points of TZ-y via the implicit function theorem, see below; we stress the implicit 
function theorem requires a nondegeneracy condition on the map. 5 

If the nondegeneracy conditions are satisfied and there is a common zero 
eigenvalue of R 1 for all 7 through x (zero eigenvalues of i? 7 correspond to unit 
eigenvalues of M~), this corresponds to the required near-zero invariant sections, 
and hence to invariant manifolds A e isomorphic and e-close to Ao- 

Let us make more precise the relation between zeroes of TL 1 near the trivial 
zero x and the zero eigenspace of i? 7 , for a fixed 7. The following lemma follows 
at once from the implicit function theorem. 6 

Lemma. Let TZ~ and i? 7 be as above; write V — T X ~E X ~ S x ~ R n - fe . Denote 
the kernel of i? 7 as K C V. Assume there is an invariant complementary space 
T. Then if K ~ R m , there is a m- dimensional manifold K C Si of zeroes for 
IZy, and T X 1C = K; for e sufficiently small, all zeroes oflZ-y within e from the 
trivial zero x lie on /C. 

Remark 2. Given any loop 7 through x, there is a loop 7' which is (nomotopic 
to the one) obtained by going k times round 7; this means that h x will include 
both M 7 and all of its powers. If fi is an eigenvalue of M 7 , there will be an 
eigenvalue \i k of My ; in particular the presence of an eigenvalue \x = e %e of unit 
modulus in the spectrum of M 7 implies that there will be maps My = (M 7 ) fe £ 
h x with an eigenvalue // = e lke . If 8/2tt is rational, there is k such that // = 1: 
thus we will have actually to require that /i ^ exp(i2irm/n) for all m,n 6 Z. 
It should be stressed that if 9/2ir is irrational there will however be k, k' such 
that \k9 - 2k'ir\ < S, i.e. \e tke — 1\<8, for any 5 > 0. This means that in 
many contexts - in particular when discussing bifurcations |18| - we should also 
require that all eigenvalues satisfy |^| 7^ 1 (see e.g. [5] for further detail). 

The maps yV( 7 leave U £ invariant; hence there are matrices M^ and M 7 ^ 
such that 

M,= (Mf -f). (I, 

We denote as i? 7 the restriction of the map i? 7 to VJ> By the above formula, 

R° := I- Mf . (2) 
In this case it is immediate to see that the kernel of i? 7 is provided by 

y_= (R*)- 1 z , (3 ) 

5 We also stress that the argument based on the implicit function theorem will give sufficient, 
but not necessary, conditions. 

e If we see Rj as a map between local sections, we should use a suitable version (infinite- 
dimensional, applying to a space of local sections) of the implicit function theorem; see e.g. 

m 
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provided the inverse (i? 7 ) 1 exists (this, of course, is the same nondegeneracy 
condition which allowed to use the implicit function theorem). The condition 
for this is just that the spectrum of M 7 ^ does not include one, i.e. that all the 
restricted characteristic multipliers (exponents) satisfy //, ^ 1 (/% ^ 0). Note 
that this must hold for all 7 (see also Remark 2 below). 

Note also that if m\ 1] = 0, we are reduced to the invariant manifold Ao 
itself, for all values of z. 

It should be stressed that in general, formula (3) will give different y for the 
same z when we consider different nonhomotopic loops 7 (a relevant exception 
is provided by the case where h x is abelian). 

Thus, at difference with the abelian case, we expect that the condition on 
spectra of monodromy matrices will not suffice to ensure the persistence of 
invariant manifolds; instead, they will have to be complemented by a "compat- 
ibility condition" ensuring that solutions to (3) for different loops coincide. 

On the other hand, if there exists an invariant manifold A', intersecting S x 
at A' x — £, the point £ is a fixed point for Af 7 for any loop 7 through x; thus, if 
x £ 7 n 7', the commutator [M 7 , My] must vanish when applied to £ S S x : the 
£,(y,z) identifying invariant manifolds must belong to Ker([/i s , h x ]) C V. 

A possible approach, employed in example II below, to discussing if solutions 
to (3) for different loops 7 are compatible is as follows. 

Subdivide Ao as the union of regions Bi, i = l,...h, each of them with 
homotopy group -K\(Bi) = Z; and denote by r\i a homotopically nontrivial loop 
in Bi , so that the homology of Bi is generated by r\i . We also write Mj for M m , 
Ri for R Vi , and so on. 

Consider then (3) for 7 = 7^: provided R® is invertible this identifies, for a 
given value of z, an invariant manifold A t over each of the B t . 

A necessary condition for the existence of an invariant manifold A e is then 
that Aj = Aj over Bi n Bj, i.e. that the invariant manifolds determined in this 
way by r\i and rjj do coincide over the intersection of the charts Bi and Bj . 

If this condition is satisfied over all the nonempty intersection regions BiDBj , 
thus determining a possibly invariant manifold A(z), we should still check this 
is invariant under the other generators of the full homology of Ao, i.e. under 
the monodromy maps for loops 7 which are not homotopic to a combination of 
the r\i for i = 1, h. 

5 The PLN theorem for involutory systems of 
vector fields 

We have now completed the geometric construction needed for the Poincare- 
Lyapounov-Nekhoroshev theorem in the general case, i.e. for involutory systems 
of vector fields Q which are not necessarily a Lie algebra, nor necessarily abelian. 

We took care to provide definitions and introduce notations such to have a 
statement which looks quite similar to the original one by Nekhoroshev EI 
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and a proof (actually, constituted to a large extent of the discussion conducted 
so far) quite similar to the one for the abelian case |17|. 

Theorem 1. Let M , Q and T> be as in section 2 above. Let Ao and U be as in 
sections 2 and 3, with Aq compact, connected and k- dimensional, and V regular 
and of dimension k in TU . Assume moreover that U is foliated into Q -invariant 
manifolds U £ ~ V of dimension n p , with k < n p < n, with Ao C Uo, and U £ 
transversal to tt^ 1 {x), so that V admits the decomposition V — ®V^' , and 
the monodromy matrices can be written in the form (1). Let rji,...,r]h be cycles 
generating the homology of Ao, and Mj = M.^ the associated monodromy maps; 

Mj = M Vj the associated monodromy matrices, with their restriction to 

the space V^°\ 

Then the following are equivalent: 

(i) There is a Q -invariant manifold A £ C U £ , isomorphic to Ao and \e\-near 
to A, for any e with \e\ < e* ; 

(ii) (a) // the loops rji and rjj intersect in some point, the commutator 
[Mi, Mj] has a nontrivial kernel inV x ~ T X T, X ; (b) The spectrum of any product 
of the matrices M^ associated to loops rjj such that M.j is not the identity does 
not include points on the unit circle S 1 £ R 2 ~ C, for all j — 1, .., h. 

Proof. It is quite clear that going around loops such that the associated full 
monodromy map reduces to the identity will have no effect on any consideration 
to follow (which justifies the specification given in point (ii)-b above), so we can 
assume all loops we consider have M. ^ L. By the Ambrose-Singer theorem, 
this excludes in particular trivial (i.e. contractible) loops. 

Let us first consider loops 7 given by rjj o ... o rjj — rjj 1 for some loop rjj. In 
this case the monodromy matrix is 

[MfY Z(j,m)\ 

/ ; 

with Z{j, m) a matrix whose explicit expression (which could be given in terms 
of M^ and M- ) is not relevant here. Thus 

R = (l-[Mf ] r -Z{j,m)\ , = (Rf } -Z(j,m)\ 

7 v oy'vo J 

The fixed points under 7 will be the kernel of this, i.e. - with the notation used 
in section 4 - will be given by 

y = [RfY 1 Z(j,m) z . 

In order for to exist for all m, we must require that eigenvalues of 

M^ are bounded away from the unit circle, so that the condition given in the 
statement is surely necessary for the existence of A E . 
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Let us now consider more general loops 7. Any loop 70 is nomotopic (suitably 
choosing the base point, see also Remark 1) to a loop 7 which can be written as 



7 = rj n o ... 

for some sequence of v € (1, n). The associated monodromy matrix is 

m, - (« J) - (; ?) 

with /? a complicate expression we do not need to write explicitly. 

As usual, we look at the matrix i? 7 = I — My, and we are concerned with 
the invertibility of Ry*\ which is now given by := I — a. Thus, in order to 
have an invariant manifold we have to require that the spectrum of a does not 
contain the unity. This is precisely (ii)-b. 

Finally, note that we need invariant manifolds associated to different loops 
going through the same point do coincide; this is precisely (ii)-a, and the 
proof is complete. <^> 



Corollary 1. Under the assumptions and with the notation of theorem 1 above, 
assume moreover that the Q -invariant manifold Ao has trivial homotopy group, 
tti(Ao) = /. Then there is e* > such that in the tubular neighborhood of Aq of 
radius e* there is a Q -invariant manifold A, foliated into Q -invariant manifolds 
A £ isomorphic to A . 

Proof. If 7Ti(A) = {e}, all loops 7 in A are contractible; hence H x reduce to 
7i°. From the Frobenius integrability of the distribution Q it follows that V is 
flat, and the Ambrose-Singer theorem |5J|^] guarantees that TL = I . (} 

Let us present some short remarks on the results obtained above. 

1. It might be appropriate to stress that in theorem 1 the condition (ii)-b 
on the spectrum of monodromy matrices is (once (ii)-a is granted) suffi- 
cient but not necessary to guarantee the existence of a continuous family of 
^-invariant manifolds, i.e. (i). To see this, just think of the case where the 
monodromy operators of all loops 7 are just My = I. 

2. We note that when A is of codimension one in Uq (e.g., if it is of 

codimension one in phase space), My are just numbers, and it is very easy 
to check the conditions given in theorem 1. Similarly, if it happens that the 
M^p commute (albeit the full M Vj may not commute) it is easy to check that 
condition. 

3. Note that if Ao is the topological product Ao = Ao x B of a topologically 
nontrivial manifold Ao and a contractible manifold B. then by lemma 3 it suf- 
fices to consider Ao (as it also follows from the role of homology groups in our 
discussion). If Ao is contractible, then our result is trivial. 

4. It should be stressed that our discussion encompasses cases where - by 
topological reasons - any vector field on A necessarily has fixed points, so that 
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the considered invariant manifold is necessarily not minimal. By a naive parallel 
with the torus case one could think that in this case the invariant manifold 
breaks down under perturbations, but actually the same topological constraints 
guarantees its persistence. In a way, only degenerations which are not enforced 
by topology are dangerous for persistence. 

5. Finally, we note that the abstract results obtained here are qualitatively 
equivalent to those holding for tori and commuting vector fields [Tl IHl fTTl ITS) : 
but now checking that the hypotheses hold in concrete situations will in general 
be quite more difficult than in the abelian case. 



6 The coordinate picture 

We have so far conducted our discussion in rather abstract terms, in order to 
make clear the geometric content of our result. However, in order to use it in 
concrete situations it is convenient to have a formulation in local coordinates as 
well; this is the aim of the present section. 

We will now introduce local coordinates in U and in Uq. Consider a local 
chart Bi C Ao with local coordinates <p\ we naturally associate to this a chart 
A{ = 7r (Bj) in U. The natural local coordinates for this will be (ip, w) where 
<p are local coordinates in Bi, and w — (w 1 , ...,w s ) are coordinates on the fiber 
it- 1 [x(ip)]. 

It is convenient to choose these as w = (y,z) where y = (y 1 , y p ~ k ) are 
coordinates in a x and z = (z , z a ) are coordinates in the parameter space 
Q. (For ease of later notation, we choose coordinates (y, z) on T, x such that 

\y%\z a \ <s.) 

The vector fields X a generating Q will be written in these coordinates as 



l dip l Ja dw k ' 

The linearized connection will be generated by the linearizations of the X a 
around A; this amounts to replacing (3 l a and /* by their linearization 7 at A: 



X° = 



0£fo,o) + o)) 



dp 1 



d 



As stressed above, we can just consider the action of the linearized holonomy 
group at x, made of M 7 for all loops 7 in Aq. Its explicit construction is quite 
standard, but we discuss it briefly below in order to fix notation. 

The nonlinear connection V provides a lift of the ordinary derivative d v along 
7, which is in general a nonlinear operator (and thus not a proper covariant 
derivative). This is written in the (ip,w) coordinates as 



7 We recall that linearization is always to be meant in transversal (to Ao) sense; this means 
linearization in the w = (y, z) coordinates, but not in the ip ones. 
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where A 1 : T —> R is in general a nonlinear function of w; as Ao is ^-invariant, 
it must be A l (ip, 0) = 0. Linearizing around 7 we have a covariant derivative 



where the matrix L, given by L* := (dA % /duP)((p, 0), is a function of ip alone. 
The evolution of j(t) = (ip(t),w(t)) along the loop coordinate t, i.e. along 
if = tpo + 2nt is described by (p — 2n, w — Lw; needless to say, the solution to 
these is 



ip(i) = ipo + 2nt ; w(t) = exp 



£[v(r)]dr 



Wq 



The function ip(t) is (as obvious by construction) periodic of period one; 
as L = L(ip), the matrix L is also t-periodic of period one, and we are thus 
considering a time-periodic vector equation for the vector variable w(t) G R' l ~ fe . 
We can then invoke Floquet theorem (see e.g. |2 111261 15U| ). 

Proposition. (Floquet theorem) The fundamental solution matrix ^(t) for 
£ = L(t)t; with L a matrix T -periodic in t can be written as ^(t) — P(t) exp[£?i] 
with P(t) a T -periodic matrix, and B a constant matrix. In particular, ^(t + 
T) = ^(t)M, with M = exp[BT]. 

The matrix M in the statement above is precisely the monodromy matrix for 
the given loop 7 C A. The eigenvalues fx±, fx n -k of M are the characteristic 
(or Floquet) multipliers; if i/j are the eigenvalues of B we have [ii — e ViT , and 
the i/i = log(/iiT) are the characteristic (or Floquet) exponents. 

The information provided by the spectrum of M is best used by passing to 
coordinates 77 via £ = P(t)n; in these the evolution equations read 77 = Br], so 
that r](t) — e Bt T]o and of course r/(T) — Mt]q. 

It is worth stressing that the unit vector field Y tangent to ^(t) in t can 
always be written in the form Y = oti(t)Xi, but in general it is not possible to 
find a 7' ~ 7 such that one can choose the oti(t) as constant, contrary to the 
abelian case. (However the loop 7 can be deformed into a piecewise smooth 
path 7' homotopic to 7 on which the cti(t) are piecewise constant.) Thus in 
the case where Q is a Lie algebra, Y does not belong to the algebra Q, but to 
the module over C°°(Y) generated by Q. In the terminology of field theory, we 
should consider the gauge algebra modelled on Q. 

With the notation introduced earlier on in this section, we will write 
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V y = d v + f l (<p,y,z)— + h a (<p,y,z)— . 



The ^-invariance of A entails that 0,0) = h a (<p, 0, 0) = 0. Linearizing 

around 7 we get 
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where the matrices F, G, H, K are functions of tp alone, and of course we have 
defined these by Fj(<p) := (df'/dy*)/^, G l a {^) := (df/dz a ) Ao , H?(<p) := 
(dh a /dyi) Ao ,K2( V ) := (dh a /dz b ) Ao . 

In other words, the evolution of the coordinates (tp, y, z) with the curvilinear 
coordinate t along the loop 7 is given by 

1 

= Fj(<p)y1 + Gi(ip)z a 
-- H?{<p)y* + K^)z b . 

The ^/-connection leaves all the level sets z — c also invariant. This implies 
in turn that h a ((p; y, z) = 0, hence M. 1 : a x — > a x . At the linear level, we have 
H?(ip) = = Kg(ip) for all loops 7, and M 7 :V ^V . 

In other words, for any 7 we can write the monodromy matrix in the form 

M 7 = [ M ? ] V? 



/ 



Hence, q of the Floquct multipliers will be trivially given by [n = 1, i = 1, q. 
The corresponding Floquet exponents are Vi = 0. We can thus, under the 
local foliation hypothesis, restrict the matrix M 7 to the (p — fc)-dimensional 

space Vo, obtaining the matrix M 7 °^ considered above; we refer to this as the 
restricted monodromy matrix for the loop 7. Its eigenvalues oti will be 
called the restricted characteristic (Floquet) multipliers, and 0i = log(«i) will 
be the restricted characteristic (Floquct) exponents. The restricted multipliers 
(exponents) do of course encode all the nontrivial information. 



7 Example I 

Let us consider V = R 3 with spherical coordinates (<p, 0,r), Q = R with coor- 
dinate A, and vector fields 

X := h(r, A) 3 V + f 2 (r, A) dg + f 3 (r, A) d r , 
Y := [a 1 (X)+a 3 (X)f 1 (r,X)]d !p + 

+ [a 3 (A) + a 3 (A)/ 2 (r, A)] dg + [a 3 (A)/ 3 (r, A)] 9 r , 

where all functions are smooth in their arguments. These commute for all values 
of r and A, [X, Y] = 8 . Their cross product is given by 

XxY = -a 2 (X)f 3 (r,X)d v +a 1 (\)f 3 (r,\)d e +a 2 (X)f 1 (r,X)-a 1 (X)f 2 (r,X)d r ; 

thus, the distribution Q spanned by X and Y is regular and two dimensional in 
Qo C Q provided the three components of this do not vanish simultaneously in 
any point of Q , i-c \\(X x Y)\\ ^ 0. 

8 Hence, if desired, we can see one of them as defining a dynamics in R 3 , and the other as 
a symmetry of this dynamics. 
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Figure 1: The double torus. The circles C and C represent the borders between 
regions I, II and III (see also fig. 2). 

E.g. for / 3 (r, A) = r — k(X) this is the case provided 

a 2 (A)/i[fc(A),A] - ai (A)/ 2 [fc(A),A] ^0 VA e Q . 

In this case there is obviously an invariant sphere S 2 of radius k(X) for all values 
of A such that fc(A) > 0. 

Suppose now that \\(X x Y)\\ ^ is satisfied on the sphere of radius ro, 
which we denote as Ao (and therefore in a tubular neighbourhood U of it in 
N = V x Q of sufficiently small radius e), and that for A = the function 
/3 satisfies /s(ro,0) = 0, (9 r /3)(r ,0) ^ 0, so that A is an isolated invariant 
manifold for the distribution Q. 

Our discussion, and in particular Corollary 2, guarantee that for |A| < e 
there is a ^-invariant manifold isomorphic to Ao (actually, for our simplyfying 
choice of the vector fields this will also be a S 2 sphere) , as the homotopy group 
tti(S 2 ) is trivial. 

It has to be noted that this can also be seen without making use of our result, 
as a simple consequence of the implicit function theorem applied to /s(r, A); 
which is not surprising as, after all, our results were also based on that theorem. 

8 Example II 

Let us consider V — R 3 , Q = R, and let A be the two-dimensional "double 
torus" (or "two- holes pretzel"), see fig.l; as well known, there is no way to 
have a nowhere zero vector field on it, so we will have more vector fields, still 
providing a two-dimensional distribution. 

Our purpose here is to build an example showing in explicit terms the validity 
of our general result in this special case. It will also be clear how to extend this 
example to the case of a pretzel with g holes. 

We decompose Ao into three (flat) regions Bi as suggested in fig.l; two 
of these are "tori with a cut", while the central one is a cylinder. This is 
also illustrated in fig. 2, where it is also shown how the three regions are glued 
together. (One could as well decompose it into two regions isomorphic to B\ and 
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Figure 2: The decomposition of the double torus into three fiat regions. The 
letters show which lines should be identified. 



-B3; we use the three-decomposition as this makes more clear the monodromy 
computation below.) With a slight abuse of notation, we will refer to the Bi 
equipped with a coordinate system as "charts" , although each of these is actually 
the union of charts, and to their union as an atlas. 

We take coordinates (xi, y^ on each of the Bi. We can take the origin of the 
coordinate systems at the center of the rectangle representing each region; we 
take Xi, yi to range from — 7r to ir for i = 1, 3, and state for ease of discussion - but 
with no loss of generality - that for the central chart x 2 G [—1, 1], 2/2 G [— n, t]. 

The Bi are immediately extended to charts (in the same sense, i.e. with the 
same abuse of notation, as above) and to an atlas on a tubular neighbourhood 
Uo C V x Q of Ao of width e; the coordinates on the chart Ai built over Bi will 
simply be {x u yi, z*; A). 

We define two vector fields on each of the Ai , given in local coordinates 9 by 

d d d d 

Xi := + tp l (x l ,y i ,z l ;X) — ; Y l := — + ipi(xi,yi, Zf, A)— . (4) 

Note that in order for Ao to be invariant under these, we must require that ipi 
and tpi vanish on A , i.e. that <pi(xi, yi, 0; 0) = 0, i>i(xi, yi, 0; 0) = 0. 

If we linearize in the "vertical coordinates" Zj and in A, with fi(Xi,yi) := 
(dtpi/dzi) Ao , gi(xi,yi) := {dipi/dzi) Ao , and similarly f l {x l ,y l ) := (d(pi/d\) Ao , 
gi(xi,yi) := {di> t /d\) Ao , we get 



A particularly simple but nontrivial choice (in which the vector field are 
chosen to be linear in z and A), corresponding to fx = f\ = ,fz = fa = 0, 

9 Note that more precisely we should add to each vector field a smoothing factor, being 1 
out of the transition regions, and smoothly sending the vector fields to zero at the border of 
their domain of definition. However, introduction of these would merely add to commutators 
some factors proportional to vector fields, so the new terms would not change the module 
structure: for ease of notation (and computation) we will just drop these. 
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/2 = q(x), f 2 = ~q(x)b 2 , gi = g 2 = 93 = 1, 9i = ~h, is the following: 

Xi := A + k,2 [q(x){zi - bi\)] A Y t := A + - 6 4 A] A . (5) 
axi 9zi dyi dzi 

Here the &i are real constants, while q(x) is a smooth function. 

Let us now consider the transition regions Ai n Aj . For the sake of concrete- 
ness - and in order to introduce a notational simplification - let us just focus 
on A\ n A2. We will write x\ = £, y\ = r/, z\ — £; x 2 — x, y 2 — y, z 2 — z. 

It is quite clear that we can just take £ — z, and that £ = £(x,y, z), t] = 
r](x, y, z) can be taken to be independent of z; so z — ( (and of course the 
parameter A) will just drop from our discussion of transition functions: these 
can be discussed in Bi n Bj . 

It is convenient to use polar coordinates (p, 1?) in B\. For a point (£,77) — 
(a;, y) € B\ n £?2, with our choice for the x range and origin, p = (£ 2 + 77 2 ) 1 / 2 is 
just p = r + (1 + x), where r is the radius of the excluded circle in B\ (recall 
< r < 1); as this function will appear often in the following, we will denote 
it by h(x) := 1 + tq + x. 

As for the angle, we can just take 1? = y. Combining these with £ = pcos(^) 
and 77 = psin(i?), we obtain at once the direct and inverse transition functions: 

£ = h( x) cos(y ) ; r) = h(x) sin(y) 

x = \/£ 2 + r? 2 - (l + r ) ; y = arctan(?7/£) . 



With the choice (5) for the vector fields, and the above, we get 



X 1 := [cosfo)] £ 
X 2 :=-§- + [q{x){z-b 2 \) 



info) 



h(£c) 



F i: =[sin(2v)]^ 
JL v — JL 

9z' ^ 2 ' 



^2 := £ + U - M 



s fo) 



dz 



+ [z- h\] 



dz ' 



Let us now consider the commutation relations. It is easy to see that the 
condition (no sum over repeated indices from now on) 

[Xi,Yi\ = a^Xi+a^Yi 

which must be satisfied in each of the Ai\[Uj^i(Ai n Aj)] for the vector fields 
to be in involution, actually imposes — — 0, i.e. = (so that 

again we can see one of X and Y as defining a dynamics, and the other as a 
symmetry of this dynamics). It is also easy to check that these are satisfied 
with our choice (5) for the vector fields. 10 

Let us now discuss the commutation relations in the transition regions: in 
each of these four fields are present, and they should be in involution. Consider, 

10 More in general, [Xj,Y;] = requires, with reference to (4'), that (dfi/dyi) = (dgi/dxi) 
and figi — gifi = (dfi/dyi) — (dgi/dxi). Defining on each chart the one-forms a; := fidxi + 
gidyi (associated to derivatives in z) and a; := fidxi + <?;dj/; (associated to derivatives in 
A), these are also rewritten as da; = and da; = a; A Note that as the Bi are not 
contractible, we are not required to have a; = d$;; actually the most interesting case will be 
the one where ai G H 1 (Bi). 
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Figure 3: The four cycles, providing a basis for the homology of the double 
torus. 

for definiteness, A\ C\A 2 . By the Frobenius condition, we must require e.g. that 
there are functions ai(x,y), /ii(x,y), smooth in A\ R A 2 , such that 

[X U X 2 ] = S := a 1 X 1 + <j 2 X 2 + ^Y, + fi 2 Y 2 ; (6) 

similar equations also hold for the other commutators. 

One can check by explicit computations (see the appendix) that with our 
choice (5) for the vector fields, eq.(6) and those for the other relevant com- 
mutators admit a well defined solution under the condition that x( x ilj) '■— 
h(x)[q(x) — sin(y)] does not vanish in B\ n B 2 , i.e. for r < x + 1 < r l5 equiva- 
lently -1 < x < -1 + 5. 

As h(x) := (1 + x + tq) > in B\ (~l B 2l we have to require that in this region 
q(x) ^ cos(y). We can e.g. require that for — 1 < x < — 1 + S the function q(x) 
satisfies |<?(a;)| > 1. This leaves ample freedom of choice for that function. 

Finally, we note that if fi,g% in (4') are not zero, then there is no manifold 
near to Ao (and homeomorphic to Ao) which is also invariant under these vector 
fields (this condition is relevant, in particular, if we are interested in bifurcations 
from invariant manifolds; see also the discussion in JS| for bifurcation from 
Poincare-Lyapounov-Nekhoroshev invariant tori) . 

In order to illustrate our result, we have to consider the monodromy matrices 
associated to cycles providing a base for the homology of A . We consider the 
cycles in A illustrated in figs. 3 and 4. 

Computation of the monodromy matrices for the cycles 771, rj 2 , 773 is immedi- 
ate: each of these lies on a single region Bi , and moreover involve only Yi . Along 
these paths, parametrized with t € [0, 27r], we have (dropping the subscript i) 
dx/dt = 0, dy/dt — 1, dz/dt — z — bX, dX/dt = 0; the solution for z(t) (with 
z(0) = zq) is given by 

z(t) = e 4 "* ^ + b (1 -e'-* ) A . 

Therefore, writing k = exp(27r), the monodromy matrix Mi (coinciding with 
the monodromy map A4i, as this is linear) and the associated Ri are given by 

_ / K (l-K)feA _ fl-K (K-l)bA 

Mi - {0 1 ) ' Rl - \ ) ■ 



16 




Figure 4: The four cycles, providing a basis for the homology of the double 
torus, here shown using the same representation as in fig. 2. 

The invariant local manifolds (under transport by V = Vo along rji) are pro- 
vided, as obvious from the explicit form of Yi, by z = biX. Needless to say, these 
local manifolds glue together into a global manifold if and only if b\ = 62 = 
b 3 = b. 

Let us now consider the path 774 and the associated monodromy matrix M4. 
The path 774 is decomposed as T4 o ... o n, see figure 4. With our choice of X%, 
only the parts n and T3 will contribute to the monodromy map; that is, z = 
for T2 and r 4 . As for n and r 3 , here x = 1, z = — 6A]. We have therefore 

to solve dz/dx = q{x)[z — bX]; writing in full generality q(x) = Q'(x), this has 
solution 

z(x) = bX + exp {Q(x) — Q(x )] (z - b\) 

for initial datum z(xq) — zq. Note this, like the equation for z{x) itself, is 
independent of y. 

It is easy to see from the above expression and simple algebra (or directly 
from the y-independence) that the monodromy map M4 : (zo,X) 1— > (z2,X) is 
the identity. (It should be stressed that this is true not only of the linearized 
map, but of the full monodromy map, whenever (p(x, y, z; X) = ip(x, y + ir,z; A) 
or does not depend on y. See, in this respect, the third remark at the end of 
section 6.) 

In conclusion, we have checked that - as actually obvious from the form 
of the vector fields - the smooth family of manifolds identified by z = bX is 
invariant under Q. 

9 Example III 

We will now provide a framework where the situation considered in example II 
is met in practice. 

Hamiltonian systems with nontrivial topology of the relevant energy mani- 
fold have been studied by a number of authors, see e.g. j2Hl; and jT^] for a recent 
contribution focusing on isochronous hamiltonian systems; isochronous systems 
on Riemann surfaces extremely robust under perturbations have been consid- 
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ered by Calogero Here we discuss simpler systems; the quantum version of 
these is studied in 

We consider a system (not necessarily hamiltonian) describing a point par- 
ticle in i? 4 , with cartesian coordinates (x, y, z, w). We also introduce "bi-polar" 
coordinates (r±, #1, r2, $2) by 

x = T\ cos$i , y = T\ sin-i^i , 
z — r-i cos $2 , w = T2 sin $2 • 

With tj = yJTj and a; = p%, y — qi, z — P2, w — (72, we would be in a 
hamiltonian framework and (I, $) be action-angle coordinates. 
In these terms, the phase space is described as 

V = T 2 x R+ . 

Let us now introduce in R a solid cone C of angle <&, with vertex in the 
origin and surface dC described in the bi-polar coordinates by 

+ tf 2 = $ 2 . (7) 

This will be a solid surface on which the point particle bounces elastically. Thus, 
the accessible phase space V a = V\C will be described in bi-polar coordinates 
as 

V a = T| X R^ , 

where T| is the torus ($1, $2) £ S 1 x S 1 without the two-dimensional disk D| 
of radius $, see (7), T| = T 2 \L>|. 

Note that each time the particle hits on the surface dC, the component po 
of its momentum in the direction orthogonal to the surface is reflected into 
— po, while other components of momentum - and a fortiori its position - are 
continuous. The motion takes then place normally until next hit on the surface, 
when again po is reflected into —po and so on. We can thus represent the motion 
as taking place on a double covering M. of V a , the two sheets of this Riemann 
surface 11 merging precisely on dC. Needless to say, this construction is just a 
manifestation of the classical Schwarz reflection principle. 

We would now like to present some remarks concerning this construction. 

• (a) First of all, note that the dynamics is singular - and actually not 
uniquely defined - on Sq :— dC; indeed, on these points the dynamics is 
defined by continuity and thus in particular it depends on the direction of 
approach (from one or the other of the two sheets of the double covering 
M). 

• (b) In view of our construction, it is entirely natural to consider a two- 
charts decomposition of Ao, each chart being isomorphic to T a and cor- 
responding to one sheet of the double covering A4; see fig. 5. This also 

llr The complex structure of this can be chosen in a number of ways, e.g. by introducing 
complex coordinates = x + iy, £2 = z + iw, or even seeing this as the product of R^j. by a 
complex surface with complex coordinate r\ = $1 + ii?2- 
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I) 



Figure 5: The decomposition of the double torus as two copies of V a , see text. 




Figure 6: The four cycles providing a basis for the homology of the double 
torus, chosen in the natural way for the framework of example III (see text) and 
represented as in figg.3 and 4. 

prompts for a different choice of the basis cycles for the homology of Ao, 
see fig. 6. 

• (c) Our construction implies the vector fields to be considered do naturally 
satisfy an antisymmetry condition with respect to Sq. This simplifies 
the computation of monodromy matrices, in particular if we adopt an 
appropriate choice for the basis cycles, each of them lying in a single chart 
of the double covering: indeed, it suffices then to compute the monodromy 
for cycles belonging to a given chart. 

It suffice now to consider a simple "unperturbed" system such as 

h = a(n,r 2 ) 
r 2 = (3(ri,r 2 ) 
i?i = wi(ri,r 2 ) 
d 2 = oj 2 (n,r 2 ) 

to get, with simple hypotheses on a,f3,ui,ui 2 - e.g. that there are pi,p 2 such 
that a(p 1 ,p 2 ) = = (3(pi,p 2 ), and Wi(pi,p 2 ) ¥= 0, w 2 (pi,p 2 ) + - that 
the system admits a double torus A in M. as an invariant manifold. This 
corresponds to the situation discussed in example II. 

If we then smoothly perturb the system, it is natural to ask if A is somehow 
preserved (upon smooth deformation) in the perturbed system. Our theorem 
allows to answer this question. 
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Finally, it should be stressed that the system considered in this example 
could be hamiltonian, and more specifically a hamiltonian perturbation of a 
hamiltonian integrable system - or also, staying within the original framework 
of Nekhoroshev's theorem pQ , a partially integrable hamiltonian system. 

In this case the unperturbed system would preserve all double tori being the 
double covering of (the part of) an invariant torus in V a , and we would be in 
the standard case of perturbation of an integrable or partially integrable system; 
however - as well known - the impact conditions would cause the system to be 
generically chaotic on the invariant double torus. 

With our construction we are able to deal with the case of perturbation of 
integrable systems with impacts, at least for what concerns preservation upon 
deformation of the invariant double tori (the construction can also be generalized 
to more complex situations). 

Note that in this case one should just compute monodromy on the standard 
invariant torus (see point (c) above); hence in this framework we just extended 
the validity of the Poincare-Lyapounov-Nekhoroshev's theorem to a class 

of systems with elastic impacts. 

Appendix. Explicit formulas for example II 

In this appendix we provide explicit solutions to equation (6) and similar ones 
for other relevant commutators in the transition region A\ PI A 2 (similar ones 
hold for A2 D A3). We recall that we have chosen the vector fields to be given 
by formula (5). 
We define 

S := (T1X1 + (T 2 X 2 + (I1Y1 + H2Y2 
(it turns out we can choose e.g. fi 2 — 0) and, for ease of writing, 

v( x ,v) : = ?0) -sin(y) ; x{x,y) :=h(x)r)(x,y) . 

We look first for solutions to [X% , Y 2 ] = S. With our conventions, a solution 
to this is provided by 

cos(y) 1 q(x) 

°"i = —, 7 > CT 2 = -, r , Mi = —, 7 , M2 = ■ 

r)(x,y) r][x,y) V{x,y) 

The equation [Xi , X 2 ] = S is satisfied with 

a l = - [h(x)q'(x)cos 2 (y) + sin(y) - q(x) sin 2 (y)] / x(x,y) , 

cr 2 = [cos(y) {h(x)q'{x) + sin(y))] / x(x,y) , 

Mi = - [sin(y) cos(y) (h(x)q'(x) + q(x))] / x(x,y) , 

M2= 0. 

The equation \Yi, X2] = S is satisfied with 

at = - [cos(y) (sin(y) (h(x)q'(x) + q{x)) - 1)] / X (x, y) , 

g 2 = - [cos 2 (y) - h(x)q' (x) sm(y)} / x(x,y) , 

Mi = [q(x)cos 2 (y) - h(x)q'(x)sm (y)] / x(x,V) , 

M2= 0. 
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Finally, the equation [Y\, Y 2 ] = S is satisfied with 

Oi = -1 , t7 2 = ^1 = Hi = . 

Note that all of these are smooth in £?i n S2 , and hence in Ai n A 2 , provided 
x(x, y) := — sin(y)] is nowhere zero in B\ n £?2- 
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